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Abstract 

Motivated by the experimental measurement of the decay rate, T, and the 
longitudinal polarization, Pl, in the Cabibbo favored decay Df — > 4>p + , we have 
studied theoretical prediction within the context of factorization approximation 
invoking several form factors models. We were able to obtain agreement with 
experiment for both T and Pl by using experimentally measured values of 
the form factors A® s<p (0), A^ s</> (0) and F D ^(0) in the semideptonic decay 
— > <pl + vi. We have also included in our calculation the effect of the final 
state interaction (fsi) by working with the partial waves amplitudes S, P and 
D. Numerical calculation shows that the decay amplitude is dominated by S 
wave, and that the polarization is sensitive to the interference between S and 
D waves. The range of the phase difference 5sd = $s — <5d accommodated by 
experimental error in Pl is large. 

PACS index: 13.25.Ft, 13.25.-k 



I. Introduction 



The branching ratio and the longitudinal polarization in Df — > <pp + have now 
been measured: 

B(D+ -> <f) P + ) = (6.7 ± 2.3)% [1] 
P L (Dj ^<pp + ) =T L /r = (0.370 ±0.035 ±0.038). § (1) 

Theoretically, Gourdin et al. [|J studied the ratio 

R h = B(Df - - 07T+) = 1.86 ± 0.26 ± °;» [4] (2) 

Within the context of the factorization scheme, which the authors of [3] adopt, 
this ratio is independent of the normalization of the form factor Ai(0). It depends 
on the ratios 

x = A^)/A^), y ee no)Mi(O), (3) 

and the q 2 dependence of the form factors. For the definitions of the form factors, see 
Bauer, Stech and Wirbel ||. No particular model for the form factors was assumed 
in [3]. Instead, Rh was studied as a function of x and y in three different scenarios 
for the q 2 dependence of the form factors. 

The result of [3] was that the (x, y) domain allowed by Rh was inconsistent with 
the measurement of x and y from the semileptonic data in Ref. f6j, and just barely 
consistent with that of Ref. [|7| . The allowed domain of x and y was also inconsis- 
tent with the theoretical prediction of 0. Ref. || also concluded that within the 
factorization scheme, the allowed range of x and y implied the following limits on the 
longitudinal polarization: 

Monopole form factors with pole-mass 2.53 GeV: 0.43 < Pi < 0.55, 

Monopole form factors with pole-mass 3.50 GeV: 0.33 < Pl < 0.55, 

Flat form factors: 0.36 < P L < 0.55. (4) 

Subsequently, authors of f| incorporated nonfactorized contributions in the decay 
matrix elements, and using the average of x and y from data [§,0, 0, showed that 
R h of (2) and 

R s£ = B(D+ -»■ <p£ + v e )/B(Dj -> 07T+) = 0.54 ± 0.10, [10] (5) 

could be understood within a scenario where the form factors have a monopole de- 
pendence as in ||. However, there had to be significant nonfactorization contribution 
to -D+ — > 07r + , though factorization need not be violated in — > <pp + . Ref. || did 
not study longitudinal polarization. 

An important point to be made is that there are three partial waves in P — > VV 
decays, S, P and D, and though the decay rate does not depend on their phases, the 
longitudinal polarization does depend on the phase difference 5s — 8d- Ref. |3| did 
not consider the effect of partial wave amplitude phases. 
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In this paper we have studied the data shown in (1) within the context of factor- 
ization invoking several form factor models (to be revealed in the next section), and 
allowing for nonzero S, P and D wave phases. This paper is organized as follows: Sec- 
tion II deals with the details and the calculation. A discussion of the results follows 
in Section III. 



II. Details of Calculation 



The decay Df — > p + <f> is Cabibbo-favored and is induced by the effective weak 
Hamiltonian given by 

H = + C 2 {uc){sd)], (6) 

where V qq i are the CKM matrix elements; Ciand C 2 are the Wilson coefficients . The 
brackets (ud) represent (V — A) color-singlet Dirac bilinears. Fierz transforming in 
color space with N c = 3 

(uc)(sd) = \(ud)(sc) (u\ a d)(s\ a c), (7) 

1 a=l 

the relevant part of the Hamiltonian can be written in the following form 

H = ^V cs V: d [ ai (ud)(sc) + C 2 8 ], (8) 

where ai = C x + ^ = 1.09 ± 0.04 @ and 8 = | £Li (u\ a d)(s\ a c). \ a are the 
Gell-mann matrices. In factorization approximation one neglects the contribution 
from the octet current part, 0§, and the matrix element of the first term is written as 
a product of two current matrix elements. It should be pointed out that there are no 
W- annihilation or W-exchange terms in D s + — > p + <j) decay. However, hair-pin graphs 
are allowed. We neglect them in what follows. The decay amplitude then takes the 
following form: 

A(D+ — = ^K s K>i (0 I sc I Dt) ( P + I ud I 0) . (9) 

Each of the current matrix elements can be expressed in terms of meson decay con- 
stants and invariant form factors. We use the following definitions: 



p + \ud\0) = m p f p e; (10) 

<<P\sc\D s > = 2 - e^e^P^V^+iiEl^mn + m^A^q 2 )- 

mo + vn^ 



£(p ' q (P» + P D ),A 2 (q 2 ) - M 2 m^ 3 (g 2 ) + 
m D + q z 

^2m^ (g 2 )}, (11) 
Q2 
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where q = Pd — P4, is the momentum transfer, f p ( for which we use 212.0 MeV ) 
is the decay constant of the p meson, e^tp) is the polarization vector of the vector 
mesons <p( p), Ai(q 2 ), (i = 1, 2, 3) and V(q 2 ) are invariant form factors defined in || . 
The decay rate is given by 

T(Dt - P + 0) = {lAx,| 2 + \A ++ \ 2 + 1^1 — | 2 } (12) 

where p is the centre of mass momentum in the final state. A 0Q , A ++ and A are 

the longitudinal and transverse helicity amplitudes given by: 

G 

A 00 (Df -> p + 4>) = -i-J=V cs V* d m p f p (m D + m^ai {aAi(mJ) - bA 2 (m 2 p )} (13) 

where the parameter a and b are defined as follow, 

l-r 2 -t 2 . k 2 

with 



2rt ' 2rt(l + r) 2 ' 

m <t> t = ^-, k 2 = (l + r 4 + t 4 -2r 2 -2t 2 -2r¥). (14) 



The other two helicity amplitudes are: 



A±±(D+ -> p + (f)) = i^VcsV^nip f p (m D +771^ j^^m 2 ) ± — 



y v H))- ( 15 ) 

The longitudinal polarization is defined by the ratio of the longitudinal decay rate to 
the total decay rate 

~r~ \a ++ \ 2 + \a^ \ 2 + \a 00 r (16) 

One can work with the helicity amplitudes or the partial wave amplitudes. We prefer 
to work with the latter as the dependence of the polarization on the partial wave 
phases is more obvious in that basis. The helicity and partial wave amplitudes are 
related by [[yj , 



A ^~Ts s+ il D ' A —T3 S± T2 P+ Te D - (17) 

The partial waves are in general complex and can be expressed in terms of their 
phases as follow 

S =| S I expz<55, P =| P I expi5p, D =\ D \ expiSo- (18) 

The decay rate is given by an incoherent sum, T oc| A ++ \ 2 + | A | 2 + | ^4oo | 2= | 

S \ 2 + \ P \ 2 + \ D \ 2 , and is independent of the phases. But the polarization does 
depend on the phase difference 5sd — 5s ~ 8d arising from the interference between 
S and D waves 
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Pl = 3^ r^-T^VTTTlT^ -■ (19) 
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S, P and D waves were calculated first by using the amplitudes ( |I3"D and (|I5D in ( |TTD 
and feeding in the phases by hand as shown in ([181). To continue with the numerical 
analysis of the decay rate T and the longitudinal polarization Pl, we have used form 
factors from six different sources: i) Bauer, Stech and Wirbel (BSWJ) ||, where 
an infinite momentum frame is used to calculate the form factors at q 2 = 0, and a 
monopole form ( pole masses are as in ) for q 2 dependence is assumed to extrapolate 
all the form factors to the desired value of q 2 ; ii) BSW/7 is a modification of BSWJ, 
where while F (q 2 ) and Ai(q 2 ) are the same as in BSWJ, a dipole q 2 dependence is 
assumed for A2(q 2 ) and V(q 2 ); iii) Altomari and Wolfenstein (AW) model fT2j, where 



the form factors are evaluated in the limit of zero recoil, and a monopole form is 
used to extrapolate to the desired value of q 2 \ iv) CDDGFN model fl3f , where the 
form factors are evaluated at q 2 = in an effective Lagrangian satisfying heavy quark 
spin-flavor symmetry in which light vector particles are introduced as gauge particles 
in a broken chiral symmetry. A monopole form is used for the q 2 dependence ( we 
mention here that we have updated this model by using more recent experimental 
results of the form factors A? K *(0), A° K * (0)&nd V DK *(0) 0, and f Ds = 241 ± 37 
MeV [14]] in calculating the weak couplings constant, of the model, at q 2 = [13| , 



which are subsequently used in evaluating the required form factors ); v) Isgur, Scora, 



Grinstein and Wise (ISGW) model where a non-relativistic quark model is used 
to calculate the form factors at zero recoil and an exponential q 2 dependence is used 
to extrapolate them to the desired q 2 ; vi) using experimental values (I| of the form 
factors at q 2 = and extrapolating them using monopole forms. 

III. Results and Discussion 

The results are summarized in Table 1 . For the entries in the last column of Table 
1 we have used the experimental values of the form factors at q 2 = 0: A^ s< ^(0) = 
0.62 ± 0.06, A? S0 (O) = 1.0 ± 0.3, 1/ Ds ^(0) = 0.9 ± 0.3 and extrapolated them with 
monopole forms. First, we note from the Table that all models, except the CDDGFN 
and the one where experimentally measured form factors are used, overestimate the 



decay rate. This fact arises from an overestimate of the form factor A±. Ref.[16] has 



noted this fact and attributes it to the imposing of chiral symmetry. Further, as Ref. 



[ IjJ has argued, more theoretical as well as experimental studies are needed for a 
better understanding of the q 2 dependence of form factors. Second, we observe that 
all the six sources of form factors allow a range for the polarization which overlaps 
with experiment with 5sd ^ 0. Note that the polarization is independent of the 
normalization of A\. It is also found that most of the final state in the decay — > 
p + cf) is in the S wave. It is also seen from the Table that the hierarchy of the partial 
wave amplitudes is: | S \ > \ P \ > \ D | . If we consider the final state to get contribution 
only from S wave the decay rates would only be reduced by (5 to 12) %, while 
the polarization would be Pl = 0.33, The hierarchy of the sizes of the partial wave 
amplitudes is in accordance with intuitive expectations based on threshold arguments. 
It is the S-wave dominance which makes an accurate determination of 5$d difficult 
(the errors in 5$d are large despite the fact that the errors in Pl are small) since the 
D-wave is an order of magnitude smaller than the S-wave. The interference term is, 
consequently, small. 
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Table 1: Decay rate and longitudinal polarization for — > P + 4>- The values of T 
must be multiplied by 10 12 s _1 . Ssd = 8s — &d is the value needed to get agreement 
with Pl data to one STD. The last column uses experimentally measured form factors. 
'Expt.FF' stands for 'Experimental form factors'. 
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CDDGFN 
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Expt.FF 
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0.32 


0.32 


0.35 




0.15 


0.37 


0.18 ±0.040 
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SD 


135 ±45 


138 ± 43 


122 ± 32 


140 ± 40 


120 ±35 


134 ±46 




s\ 

p \ 


4.3 


3.7 


3.8 




2.8 


5.5 


4.7 




s\ 

D\ 


11.9 


13.5 


7.4 




8.2 


8.6 


16.5 


Experimental values of T and Pl 


r 


= 0.14 ± 0.05 0,P L = 0.370 ± 0.052[| 
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